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1>EVEL0I'MI NT OF \ nmiE-DlMI NSIONAL \A\ IFR-STOKFS COl'F 
ON i.PC STAR* 100 v\WllTFR 


Bv 

\>*>r N. VAtait^ 

INTROPIXTION 

The Import 4ncc of <i itenoral purpose, three*illmeii«ion4l \.iv»er*StoKe« solver 
cannot he overemphasl :ed due to the Inadequdcv of houndarv-laver evpuitlons and 
other approximate equation sets to solve complex flv>»« phenomena such as those 
encountered In the vlcinltv of a i«lng»bodv junction, hlille no coie appar»*ntl\ 
exists that can treat every flow situation, especially for turbulent flows, 
significant progress has recently been made bv several Imesttgators irefs. I 
to 4) using biHindarv-fltted coordinate svstems (ref. .A) to treat comple . 
geometry. Once the governing equations are cast In the bi*dv-fltted coordinate 
system, the approach Is valid for different conflg»irat Ions and requires onlv 
the appropriate metric coefficients of t ransformat Ion from a suitable code, 
such as that described in reference A. Pue to the aval labt I ity of a special- 
purpose computer U*IH' SrVR-lOO) at NASA laRi”, which Is capable of achieving 
substantial speedup compared to conventional computers (refs. S to for 
explicit type algorithms, an effort was initiated to develop a three-dimensional 
code In bodv-fltted coordinates using Macl’ormack ' s algorithir iref. '1. The 
code Is structured to be compatible with anv general conf Igurat Ion, provided 
that the metric coefficients for the t ransformat I .'n are available. The 
governing equations are developed in primitive variables lo , u. v. w, p and T) 
in order to facilitate tiie Incorporat Ion of physical boundary conditions and 
turbulence-closure nw’dels. Mact'ormacK's two-step, unsplit, time-marching 
algorithm ivef. is used to solve the unsteady Navier-Stokes equations 
until steady-state solution Is achieved. 


' Research Associate, Old Pomlnion University Research Foundation. Norfolk. 
Virginia. .M.APS. 
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GOVERNING EQUATIONS 

The three- .Inenslonal, compressible Navler-Stokes equations In Cartesian 
coordinates are as follows: 


where 
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The street and flux terms used in the above equations arv given by 



(5) 


(4) 


Note that Stoke' s hypothesis (X ■ - j u) is used in the stress terms. The 
internal energy e is given by 


e ■ C T ♦ -^(u^ ♦ 
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Th* p*rf*ct g«f law «n«l Sutherland's vltcatltv equations are used to close the 
systesi 
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where and u< ar« constants. 

The above equations are in dinensional for* and are valid only In the 
Cartesian coordinate system lx.y,:). These equations are then transformed to 
a new coordinate system iC.n.C) using Implicit dl fferent lat Ion. Finally, the 
equations are nondimensional 1 :ed with appropriate reference quantities, 
details of these transformat ion and final f«'Tiis of the equations are presented 
In the \ppendlx. 


NINFRICAL SaiEME 

The CDC STAR-100 Is a special tvpe of computer which can do vector 
operations at a verv fast speed. The word vector Is used here to describe a 
set of up to oSOOO data elements stored In consecutlxe storage locations in 
memorv. This pai’tlcular feature makes explicit methods verv attractive for 
the STVR-lOO computer. MacCormack iref. *) developed a two-step marching 
algorithm which has been used on a variety of problems .ind Is especially suited 
to obtaining steady-state solution of the Navier-.*»tokes equations. To explain 
the algorithm, we recall (see Appendix) that governing equations can be written 
In the vector form as 
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where U, F, C, and M are vectors of length 5 and are given bv equation 
IV-ll) in the Appendix. Let superscript n denote th» variables .if a given 
time t. The variables are .idv.inced to a new time t * dt through predictor 
and corrector steps that are given as: 


Prtdictor tfp . 


"fl.k • "?j.k • • "I.j.k)- • s.).i) 

•(!?) • ''J.j.k) ■ ‘“(“?.j-i.k ■ ‘■'u.k) 

k-»;,),k)] 

■(«)[*>>(''l.),k*l • ’’l.J.k) ' *2>(°l,J,k.l • °l.j.k) 

• “»("?.J,k.. ■ "J.j.k)] <’> 

The first derivative terms used In flux vectors F, G, H are evaluated 
using two'polnt backward differences. 

Corrector step . 

“l,!,k • J [“l.J.k * “uj.k] 

•(!l) I [*“('’M.k • ''l^.j.k) • 

• »»("M,k • Cj.k)] 

•(!^)T[«.:^rI,k • C,.k) • - =fl-,,k) 

"(sc) : [®>’(''i.i.k ■ ''iil.k-l) * * 2 >(“lJ,k ■ °l,j.k-l) 

♦ “»(«M.k ■ "M.k.l)] 


( 10 ) 


where the first derivative terns used In flux vectors F, M are evaluated 

using two-point forward differences. 
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Successive applicetlon of predictor and corrector step equations advances 
the solution by time At. 


Damping Terms Used In Navler«Stokes Equations 

Fourth order damping has been used in the Navler«Stokes equations to capture 
discontinuities (e.g. shock waves) in the solution. The particular form of 
damping chosen Is from Tannehlll et al. (ref. 8) and Holst (ref. 9), where a 
source term Is added In both predictor and corrector steps of the solution. 

The main reason for selecting this type of damping is that it is easily put in 
vector form on the STAR*100 computer. The form of damping used Is given below. 

Predictor step . 


,r . I Pt.l.l.k • • Pl .t.k • Pl-M.kl 

i.Jik if n ♦ ? n'' ♦ n" \ 

PiM.J.k ^ Pi.J.k Pi-lJ.kj 

* 1 u? , , L • 2 U? , . * l/J , . . 1 

li*l»j*k iij.k i'ltjt k J 

^ (^l,j*l,k * ^ PiJ.k * Pi-l.j.k) 
“ ['^t.J^l.k ' ^ ^ij.k * ^ij-l.k] 

* (’’i.j.k^l ^ ^ Pi.J.k * Pi.J.k-l) 

'' ['^i.j.k*l ’ “ ^.J.k ^ ^i.j.k-l] 


Corrector step . 
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(concl 'd) 


Cj are the damping coefficients generally taken to be 0.1; 


Here C , C 

X y. 

however, the specific value is a function of the problem being solved. Damping 

can be completely removed by equating C , C , C, to zero. 

X y 2 


Arrangement of Data 

Becuase of its special architecture, this section is devoted to the 
arrangement of data in the St’AR-100 computer. The computer uses the "virtual 
memory" addressing system, which allows a code to have a data base larger than 
the central memorv of the computer. However, when such a large code executes, 
the system can only accommodate S24,2S8 words of memory (distributed on 7 
magnetic disks, referred to as pages); the rest of the required information 
is stored or. peripheral disks (pages). When the code refers to some infor- 
mation that is not resident in the central memory, one of the pages (disks) 
from central memory becomes a peripheral unit to accommodate the required infor- 
mation in cores, '.his process is called "paging" and is a slow process (taking 
about 1/3 sec) due to hardware limitations of the system. Due to this feature 
of the system, the data in the code is arranged so that maximum use is made of 
this information residing in the memory at all times. Lambiotte (ref. b) 
addressed this question in an earlier paper; following the procedure suggested 
in that work, we interleave all the global variables into one variable and arrange 
the data plane-by-plane from inflow plane to the downstream end. To give an 
example, let A be the interleaved global variable A(IMAX, IVAR, JMAX, K.M/\X) 
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^h^r• IMW, JtuI KM*.\ .u'f th«» mimht»r of grul polnH 1» v, <iiul 

I Jir#otlon% rmpeot 1 veb .tiul l\ \H rofcrt to lUffcronl for 

o\ampl<«, hv 4«*lgntiig l\ \K • I for tionvU^, \ilMA\, I, .IHW, kMVM ^>ulUI rrfor 
to thr il#iuttv, thin notation la nuAtlr nrnuinttc for rcAilthtIt v i».g. one IKIh> 
for driultv, ^tc.) itnJ 4ppll<>il to 4lt glohitl v.irl4hl«>« it two ttm^ l^\#ln, thi* 

4>iviint«ge of «uoh «n 4|«|«roiioh In thiit nt .1 glvon I location, .ill the r<>«)ul 
vnrlnbUn «r» grvMipi*»l togrthi'r pornlttlng v#otor oprrnt l»*nn of nUi* ,rMA\*ikM\\ 
nt #ni'h plnnv. to mlnlmlt^ thr pnglng further, i»i» perform thr prr»IU*tor 

ntep from I • 2 to IMA\ • I .iml rorri»ctor step from I • IMAX . I to no thnt 
tionr the tloKiint renrn houtuUrv (1 ■ IMAX) ko mnXe hrni tine of the .ivnlhiMe klnt.n, 
thin iipprt'iU'h hnn hreii veri niu'oennful , .iito immorloiil experiment 4t lot* i»lth the 
Ov'mptiter ooile iiulUMten th;it the reniilting pnge fniiltn were molnlv «lue to the 
vllfferenoe between vhitn bnne .mil oentr.il memorx’ iiul ootiUI be well pieilloteil 
before the .ictuiil eneoutton of oo»le, Tor 4 Al».A|»3l grid, the il.it.-i bane 
Inoliiiilng coding wnn nboiit II pngen long, nnd It remilted In .in 4\er.ige of I 
pnge fnultn ill - • 41 per Iterntlon cvcle, 

KISUl.tS \M' PISniSSION 

the remit n prenentod In thin report .ire prelimln.irx In mitiire, but 

neverlhelenn re|."enc:tt .tn Impm-t.int vl.inn of problemn th.it could be tmed .in 

tent C4nen for future Invent Igntorn involved with three dlmennlon.il \4V ler 

Sioken problemn, Vhe compiit .it loiifil domiiin npecl.illjed to .in 4vl.il corner m 

nhown In figure 1. Vhe Inflow boiindnrv in ,it t • 0 .ind outflow boiindarx .it 

■ • t . Vhe outer edge In formed bv n • n -md 5 • boiindorlen, 

m;iv mnv mnv 

wherein the w.ill boinulirv in formed bv n • 0 ind t • 0. Vhe rent oi the 
nei.* In devoted to virloiin lent problemn ind their boiindirv conditlonn. 

i‘ine l! flit riite In .’'upernoniC free Streim 

Vo tent conditlonn chonen ire • .A.O, Ke.^ • i',.-.l ' lii' per foot; 

V • liid»"*M ind Vw • O.ls.do'H. Vhln tent cme In chonen to v.ilidite the 
t 

computer code nince iccunte boiindirv l.iver renultn ire ivulible lor 
compirinon. Vhe flit plite In locited it the n • i' l.ice of the compiitit lon.il 

box of figure I Vhe following boiindirv conditlonn ire ipplled 
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M ( • 0 All th<> l'|o«t v«rUMtf« Arr At th^lr 
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fr<><**«tr»AM vAlur, 


Ha I I _ C iMuii t ton. « \t n • 0 1^ • 

u • V • • i' tno Allp. n»» I It 1 1 1 **n > 

r • rw 



\» t • 0 i.i • n, cont ImiAt l»*n coiulltloit I* u»#%l on All thr flot. tArlAblp*. 1 .#. 
\A|U(» of A VATIaM# At .1 • I l> *«*t Ai^IAI tO lt« \aIiIO At .1 • 

OutfV otliio oon>lltlon.» \t n • m Attil ; • . contlmiAtlon oon>Utton 

t» on All vArlAMcA. 

i>utflot* ooiulltton,. Vt ; • t . i*ll »>'«’ \ av 1 aM«*» .irn ol»t.«lno»l n»\n>j 

— — - . — » niitx 

irrv'th or»l**r #xtrAj»ol At Ion, »*Mch I* Alont to tin* oontInuAtlon coiulltlon. 

In fliiur# i, A\Ul volooitv r«>«ultA Ar«* *lto*.n for thr«*r Jlffrrrnt oo*I«*a. 
SUIMKMV r«*f«*rA to M.ui'ormAoK • a Ha|M»I Soltor ir^f. liM, a oo^Io r«*loA#«»»l l*» 
th«* \\S\ VmoA KrAfAroh Tontor. Othor vUtA on tUuri* I a from the un|nil*l i Aho*l 

roAiiltA of I'AVp II. Hiuh of N\:^\ IaKi'. H»*una.»rt . I ax or Aoliitlon tA aIao Alu*x.n 
In tltlA flituro, Vho oAUnil.ttion* horn >*rvp »lono mtli ^I'O'JI mpAlt, thnx onix 
.•| no.loA wotv iiAO.I noimAl to tiu* plAtc aiuI *1 n.slo' .tionj: tho lonuttiultnAl 
klirpv'tlon, VliP itormAl moAlt ia Misitlx Atmoltr.l to jjlxo jtoo,! rpAoliit (on (n the 
xlAOouA lAxor. I r\xm thlA ftjturp, It (a nototl thAt aII tltror oo»Ipa bAAlCAlIt 
jjtxp AlmilAr rpAultA. HxIa toAt oaap Atmt'lx aaa trlp»l to ijtxc Aomr v*rp>l(Mlltx 
to the nox, covIp, 

i'aap JI; SuporAOtUo ri«*x« \lon >5 .tn Vx(aI I'ornor 

rhr noxt tPAt OAAp cltOAon for tlt(A Atiulv (a of AupprAonlo floA aIohiS aii 
axIaI oornor, Vito toAt ooiulttionA .iro l.lonttOAl to thoAo of iMap I, *,nr 

rpAAon for Aplp,‘t(ttji tItiA oaap haa tlxAt I'r. f.M, Uunis ot \\S\ \moA Ixa^ AlroA.lx 
rxin th(A OAAo Aiul it aaa I«p|(oxp»I tItAt A oompAi'lAon A(tlt An (inlpi'oiulont Ovslo 
aoiiKI a.M to tlto vM-r.lU'llltx of titr jMOAput ot'.lo. 

Vito I'ouitvlArx oo*.oit*onA for tito axIaI vornrr ai'o aa tollonAt 

Inflow oonJltton. At t • 0 All tho fl >w xaimaMoa tro aaauiuo.I At tho(r 
froo-*AtvoAm vAltio, 


II 


ion . » \t n • 0 iuul ; • I't 


?r 


u • V • >« • 0 (no *ltp, no tnicct ion) 

r • u 


i' 4nii 


n 


«»ila » oonditlon .’ 


««\ 


AlUl 




I'ont Imut Ion oo'iJltlon t« vised on All VArlAPles At 


vXitflovt condition.- \t c • i . <»ll tX«* vArl.iMes Are obtAlned bv extrA- 

polAt Ion, 

It shxMild be noted th.it, in .iddltlon to the lnt‘l.*»* plAive, there Are 2 nu're 
plAiies: I • .* And I • In the free streAm, Pie corner stArts At I • 4; there- 
fore, At I • And I • AppropriAtv svmmetrx And Aiitisvmmetrx conditions Are 
imposed At n • 0 And ; • v>, K .'I'.^l'.M mesh Is used for All these comput At Ions , 

The results Are shoMii In fljjures .1 to .*«. Ml the compnrlsons Are done At 
X • v'.i>55S5', In fUures \ to e. axI.iI xeKvltx. crossfl.'w velocity, densitv 
And pressure profiles Aloti|t the dlAjionAl At x • v'.vtdS.so* .ire shov»n, Ixcept 
f,'r the cri'ssfli'w velocltx, which Is );enerAll\ the m,'st sensitive ,pi.intltv, 
compArlsons Are in reASoiiAble A){reement , \,*te Alse thAt in Munij’s scheme, 

difference equAtlons .ire written In phvsicAl sp.ice inonuniform mesIP . whereAs 
In the present Cv'de the difference e»piAt Ions Are written in c.'mput At lon.il spAce 
lunlform meshl. Such different .ipproAches h.ive different trunCAtion erriM's ,ind 
could ACCv'unt for the sll);htlv different results observed here. 

In flijures ' and velocity And densitv profiles .ire c.'mp.ired At 

X • i’,i'9SSe* And At two fixed VAlues of N • I'.i'le^* Is In the corner 

I AVer whereAS \ ^ Is f.ir awav from the corner where we expect to see a 

mA X 

boundAi'v - 1 Av er tvpe beliAVlor. I'ompArlsons sh.’w a jjood Ajjreement between the 
two schemes, vhi fl>sure ■* boundArx • I a\ er s.»lution Is Also shown. aiuI it A»;rees 
verv well with the \av ler-St oKes solutions .it > • i»s expected. I'hls set 

of results thus indlCAtes th.it the present three dimenslon.il N.i\ t er-St oKes 
code IS b.isicAllv sound aiuI woul.l >{t\e correct results f,>r de)jen''r.it e c.ises 
ve.>:, fl.it plAte And Axi.il corneiM . 



Case 3* Subsonic Flow in an Axial Corner at M • 0.9S. 

m 

The test case chosen here is of the subsonic flow in an axial corner at 
• 0.95, for which asymptotic solutions in a corner layer have been presented 
by Weinburg and Rubin (ref. 11). 

Boundary conditions .- The only difference in this case compared to the 
supersonic case is ir the outflow boundary condition at C • £ . In this 

case, pressure at the downstream end is specified at its free-stream value. 
Velocities .and temperature are obtained by extrapolation, and density is then 
computed using equation of state (A-19) in the Appendix. The wall temperature 
is taken to be equal to free-stream static temperature. 

The axial velocity and crossflow velocities along the diagonal at x • i).n' 
arc shown in figures 9 and 10 in corner- layer coordinates (ref. 11). The 
Navler-Stokes code was run at Reynolds numbers, R • 10S38 and R ■ 4.152 to 
see the effect of Reynolds number on the solution. From figures 9 and 10 it 
is clear that as the Reynolds number Increases, Navier-Stokes solution moves 
towards the as>Tnptotic solution, which in this case is the solution to corner- 
layer equations for infinitely large R^,: however it is still quite far from 
the as>'mptotic solution. On the basis of experience gained in similar studies 
with interacting boundary- layer solutions (refs. 12, 13), it is conjectured 
that a value of R ■ O(IO^) may be required to reproduce the as>Tnptotic 
solution. However, high Reynolds number flows are verv time consuming as one 
can see from the stability criterion, the Courant -Friedrichs-Lewy (CFL) 

Condition; 


At < 


1 




♦ 




113 ) 


where a is the local speed of sound. 

Note that as viscous layer becomes thin at higher values of R^^, smaller 
step sites (Ay, A:) are needed to resolve the viscous layer; hence a smaller 
value of allowable time step At results, thus taking longer time to converge. 
Because of insufficient time, higher Revnolds number cases were not tried; 
however, present results should provide a good reference point for future 
investigators who miglu be Interested in resolving the details of corner 
layer. 
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In fluure 11, axial velocity profile* at z • z^^^ are compaied with the 
tH>uniiarv«laver solution. It is clear frosi this fiijure that at R • 4^1S^, the 
Navier-Stoke* solution recover* the hounJary-laver solution. Use of the same 
set of eipiations throujjhout the computat lonal domain (corner- 1 aver boundarv- 
layer inviscid regions) and recovery of boundary- 1 aver and Invlsctd solutions in ' i 

their respective doisains assure that the comer- layer solution so generated is 

the correct solution. I 

I 

Case 4: Supersonic Flow in an Axial Corner at • l.S ^ 

The last teat case presented is that of supersonic flow in an axial corner | 

at • l.S, for which as>Tsptotic corner-laver solutions are available in 
reference 11. The wall temperature is set ei\ual to the free-stream static 

I 

temperature. IWnindarv conditions are Identical to those of test case J. j 

The axial and crossflow velocities along the diagonal at x • 0.4* are shown » 

in figures 1. and 13. The corner-laver solutions are also plotted in these ] 

figures. It is noted that although the \avler-Stokes solution comes closer to I 

the as>-mptotic solution as Reynolds number is increased, the shift i* verv | 

slow compared to the subsonic case. Rased on this limited data base, it seems 

very unllkelv that Navler-Stokes solution will merge with asymptotic solution 

at finite values of Reynolds number. As in the previous case. It Is seen from 

figure 14 that at : ■ : the Navler-Stokes solution has recoxered the 

max 

boundary- layer solution except near the edge of the boundarx laver, where we 
do not expect good agreement due to leading-edge shock in the Navler-Stokes 
solutions. 

It is felt here that the corner-laxer solutions of kelnberg and Rutin 
(ref. 11) do not appropriately model the flow in the supersonic axial comer 
since no mechanism exists that can model the complex flow pattern that develops 
in the corner laver by interacting shocks. Figure 15 shows a constant density 
contour plot at x • 0.4’ obtained from the Navler-Stokes solution. The leading 
edge shocks and the resulting corner shock are identified on this figure. Sucli 
a shock system will have considerable Influence on the viscous flow In the 
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corner. 


CONCLUSIONS AND Rt!C(»Ml:NDATIONS 


A computer code for solving three-Jlmensional Navier*StoLes equations has 
been developed for the CDC STAR* 100 vector*processing computer. The governing 
equations have been transformed through implicit differentiation; consequently 
the code can be applied to any general configuration provided the appropriate 
metric coefficients of transformation are available. The code has been tested 
and validated for several flows through comparison with known solutions. 

Results for the subsonic axial comer, while encouraging, are not complete 
and require additional study. It is recommended that studies he made at higher 
Remolds numbers in order to reach definite conclusions concerning the agreement 
between the corner-layer solution and Navier-Stokes solution. 

For the supersonic axial corner, the comparison is not that good. It is 
conjectured that lack of modeling of complex, leading-edge shock interactions 
in the corner-layer equations is responsible for the obser\’ed differences. It 
is recommended to validate the solutions of Navier-Stokes equations for super- 
sonic axial comer with experimental result) when these become available. 
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Figure 4. Crossflow veloc 
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;y along diagonal (M^ ■ 3.0). 
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Floure 9. Ax1.i1 vf'locUy .llono d1.ioon.i1 (M • 0.95). 
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APPENDIX 


TRANSF0R.HXT10N OF an’FRNINll Ev:UATIONS 


The tet of govominn equations (eqs. I to are written in the Cartesian 
coordinate system. These are now transformed to a new set of independent 
variables ((.n.C) that form the computational space by using Implicit differen* 
tiation. The new set of independent variables is completely general and forms 
a paral lelepiped* shaped computational domain. 

Let f be a function to be differentiated. 

Since 

C • C(x.y.:) 

n ■ nlx.y,:) 

C • c(x,y,:) (A-l) 


one can write 


f. ■ 

x.f 

♦ V, f * 
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♦ V f ♦ 
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n X 
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♦ V f ♦ 
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(V-5) 
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where 


• W’ * %<>’;={ ■ >'c=t> • ''«'>'{> ■ N=c’ 


(A.6) 


1< the Jacobian of transformation. Rearranging the terms, the first deriva- 
tives of f can be rewritten as 


where 


it 


}7 - 


(A-7) 


3f 


♦ 822 ^,, ♦ 82 jf^ 


(A-8) 


If ■ Sjifj ♦ ♦ 8jjf^ 


(A-9) 


^12 • (y^ij 
BJ3 • (Vf-, 
«2l • 

822 ■ 

B 23 ■ (x,:^ 

031 • cx^y^ 

032 • (x^y^ 

033 • (x^y. 


y^:,)/J 

X^z^)/J 

x^:^)/J 

x,y^)/J 

x^y^)/J 


(A-IO) 


Using equations (A-7) to (A-9), equation (1) can be transformed to 


3u 


♦ BllF, ♦ 62111^ * - 31 ^<r * ^12^^,, * 022^^ B33H, 


♦ 013^^ ♦ 023^^ ♦ 03 3^U ■ 


(A- 11) 
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Th« Shear terms In F, G. H (eqs. 5 to 4) must be transformed to the C. 
n. C system also. These are given below: 


T • 

XX 


T • 2u 

yy 




T ■ T 

xy yx 


* 8i 2U^ ♦ gljU^ - 

.•i] 

r <»11V5 ♦ ai2v, 

L*8:i“r ♦ 822 “p 


BjjV^ ♦ ♦ fl2jV 


Bjiw ♦ 6,2W ♦ B,,w 


BijV 

82 jU 


:] 

. . . „r "’'“t * *”“n ■ 

■* ^11“^ ♦ ®12“^ ♦ 

- T, . u[ * '”“0 

” L* ♦ «J32V^ ♦ ^3ivj 



(A-i:) 


A • Bn“5 * ®i2‘*n * 

♦ 82 JV^ ♦ ^22^^ ♦ 823V^ 

* ®3I**f * '*32“„ * ®33“, (A-13) 

^ n >. 


Nondimensionalization Scheme 

The equations presented so far are dimensional. However, it is generally 
more convenient to work with a nondimensional form of the equations, especially 
since important parameters such as the Reynolds number explicitly show up in 
the nondimensional form. The dependent and independent variables are nondimen- 
sionali:ed according to the following scheme: 
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lA-14) 


where * euperecripf refet* to 3 Jimeneional quantity. Xoc that thl» »uper- 
script has not been used in the equations presented so far in the discussion to 
avoid the cumbersome notation. However, it has to be made clear that the 
equations up to equation (A>13) in the Appendix and up to equation (7) in the 
main text are in the dimensionwl form. The remainder of the equations are in 
nondiswnsional form. Letting 
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0 u L 
r r r 


(A-I51 


be the Reynolds nurabec, the nondimens iona I forms of equations (2a) to (2d) are 
as follows: 
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where 
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Also. »e have 


e • 03 ! 




(A. IS) 




3T 


(A- 19) 


To fiinmariio, the govorning equations in nondimenslonal form are the some 
as equation (A-1) where the vectors U. F, G, and H are given by equations 
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(A«l 6 a) to (A*1 (hI). The shear stress terms In nonJlmenslonal fora are identical 
to equations (A*l«) to (A>13), provided a nondiaensional fora of variables is 
usvd in the definition of these quantities. 

Tho t>'pes of nondiaensional schemes are optional in the computer code: 

(1) The first scheme is based on flat plate variables and uses t.hc 
following reference quantities 

u«u, P*P, T»T 

**r " ^r " * * “r * 


This is most suitable for flows that have straight flat rejtions at 
lead inn edne. 

(2) The second scheme is based on Van PvKe's variables. Here, the 
following reference quantities are used 


u ■ u . 

r • 





This set of variables is most suitable for blunt body flows with 
bow shock ahe.!!! of the body. Any other nondimensionalltation scheme 
con be easily incorporated into the code by suitably modifying the 
constants Ci, Cj. Cj, C 4 defined by equation (A-P^. 
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